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Abstract. We give an overview of recent developments in Sturm-Liouville 
theory concerning operators of transmutation (transformation) and spec- 
tral parameter power series (SPPS). The possibility to write down the 
dispersion (characteristic) equations corresponding to a variety of spec- 
tral problems related to Sturm-Liouville equations in an analytic form 
is an attractive feature of the SPPS method. It is based on a computa- 
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—j. ' tion of certain systems of recursive integrals. Considered as families of 

functions these systems are complete in the Z/2-space and result to be 
the images of the nonnegative integer powers of the independent vari- 
able under the action of a corresponding transmutation operator. This 
recently revealed property of the Delsarte transmutations opens the way 
to apply the transmutation operator even when its integral kernel is un- 
known and gives the possibility to obtain further interesting properties 

Ky{ • concerning the Darboux transformed Schrodinger operators. 

We introduce the systems of recursive integrals and the SPPS ap- 
proach, explain some of its applications to spectral problems with nu- 
merical illustrations, give the definition and basic properties of trans- 
mutation operators, introduce a parametrized family of transmutation 
operators, study their mapping properties and construct the transmu- 
tation operators for Darboux transformed Schrodinger operators. 
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1. Introduction 

Transmutation operators also called operators of transformation are a widely 
used tool in the theory of linear differential equations (see, e.g., |3], [ID]) @H]> 
[5Tj . [531 and the recent review [HT]). It is well known that under certain 
quite general conditions the transmutation operator transmuting the oper- 
ator A = —-j^z + q{x) into B = —4-r is a Volterra integral operator with 
good properties. Its integral kernel can be obtained as a solution of a certain 
Goursat problem for the Klein-Gordon equation with a variable coefficient. 
In particular, the elementary solutions of the equation Bv = Xv are trans- 
formed into the solutions of the equation Au — \u. If the integral kernel of 
the transmutation operator is unknown, and usually this is the case, there is 
no way to apply it to an arbitrary smooth function. This obstacle strongly 
restricts the application of the transmutation operators confining it to purely 
theoretical purposes. 

Recently, in [3] a relation of the transmutation operators with another 
fundamental object of the Sturm-Liouville theory was revealed. Sometimes 
this object is called the L-basis 23] where L refers to a corresponding linear 
ordinary differential operator. The L-basis is an infinite family of functions 
{fk}kL such that Ltpk = for k — 0, 1, Lipk — k(k~ l)ipk-2, for fe = 2,3, .. . 
and all cfk satisfy certain prescribed initial conditions. In [IT], [55], [55] it 
was shown that the L-basis naturally arises in a representation of the so- 
lutions of the Sturm-Liouville equation in terms of powers of the spectral 
parameter. The approach based on such representation is called the spectral 
parameter power series (SPPS) method. The functions ifk which constitute 
the L-basis appear as the expansion coefficients in the SPPS. In [4jJ, [52] and 
|45] convenient representations for their practical computation were proposed 
which converted the SPPS method into an efficient and highly competitive 
technique for solving a variety of spectral and scattering problems related 
to Sturm-Liouville equations (see [T2], [13], [37], [39], [45], [47]). The above 
mentioned relation between the transmutation operators and the functions 
ipk called in the present paper the recursive integrals consists in the fact es- 
tablished in [5] that for every system {y>fc}fc = n there exists a transmutation 
operator T such that T[x h ] — ifk, i.e., the functions ifk are the images of 
the usual powers of the independent variable. Moreover, it was shown how 
this operator can be constructed and how it is related to the "canonical" 
transformation operator considered, e.g., in [STJ Chapter 1]. This result to- 
gether with the practical formulas for calculating the functions ifik makes it 
possible to apply the transmutation technique even when the integral kernel 
of the operator is unknown. Indeed, now it is easy to apply the transmutation 
operator to any function approximated by a polynomial. 

Deeper understanding of the mapping properties of the transmutation 
operators led us in [46] to the explicit construction of the transmutation 
operator for a Darboux transformed Schrodinger operator by a known trans- 
mutation operator for the original Schrodinger operator as well as to several 
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interesting relations between the two transmutation operators. These rela- 
tions also allowed us to prove the main theorem on the transmutation oper- 
ators under a weaker condition than it was known before (not requiring the 
continuous differentiablity of the potential in the Schrodinger operator). 

In the present paper we overview the recent results related to the SPPS 
approach explaining and illustrating its main advantage, the possibility to 
write down in an analytic form the characteristic equation of the spectral 
problem. This equation can be approximated in different ways, and its solu- 
tions give us the eigenvalues of the problem. In other words the eigenvalue 
problem reduces to computation of zeros of a certain complex analytic func- 
tion given by its Taylor series whose coefficients are obtained as simple linear 
combinations of the values of the functions ipk at a given point. We discuss 
different applications of the SPPS method and give the results of some com- 
parative numerical calculations. 

Following 9^ and 146] we introduce a parametrized family of transmuta- 
tion operators and study their mapping properties,we give an explicit repre- 
sentation for the kernel of the transmutation operator corresponding to the 
Darboux transformed potential in terms of the transmutation kernel for its 
superpartner (Theorem I6.2J) . Moreover, this result leads to interesting com- 
mutation relations between the two transmutation operators (Corollary I6.6[) 
which in their turn allow us to obtain a transmutation operator for the one- 
dimensional Dirac system with a scalar potential as well as to prove the main 
property of the transmutation operator under less restrictive conditions than 
it has been proved until now. We give several examples of explicitly con- 
structed kernels of transmutation operators. It is worth mentioning that in 
the literature there are very few explicit examples and even in the case when 
q is a constant such kernel was presented recently in [9 . The results discussed 
in the present paper allow us to enlarge considerably the list of available ex- 
amples and give a relatively simple tool for constructing Darboux related 
sequences of the transmutation kernels. 



2. Recursive integrals: a question on the completeness 

Let / £ C 2 (a,b) D C l [a : b] be a complex valued function and f(x) ^ for 
any x £ [a, b]. The interval (a, b) is assumed being finite. Let us consider the 
following functions 

X^(x) = 1, xW(s) = n f X^~ l \ S ) (/ 2 (s)) ( ~ ir da, 



Xq 



x £[a,b], n = l,2,.... (2.1) 



We pose the following questions. Is the family of functions {X™} 
complete let us say in Li(a, b)1 What about the completeness of {X( 2n ^} 



or r x (2«+i)}°° ? 

I J n— 
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The following example shows that both questions are meaningful and 
natural. 

Example 2.f . Let / = f, a = 0, b = 1. Then it is easy to see that choosing 
x„ = 0we have Jf(°)(«) = 1, X^(x) = x, X^(x) = x 2 , X^(x) =x 3 ,.. .. 
Thus, the family of functions {X^™^} is complete in £2(0, 1). Moreover, 
both {X( 2 ")}^° =0 and {X (2n+1) }^ =0 are complete in L 2 (0,1) as well. 

If instead of a = we choose a = — 1 then {X^} „ is still complete 
in L 2 (-l, 1) but neither {X^}^ Q nor {X( 2n+1 )}~ ?~ ' 

Together with the family of functions {X^™)} we consider also an- 
other similarly defined family of functions {X^} „, 

X(0) = l, X^(x) = n f X^-^is) (/ 2 (s)) ( " ir " ds, 

Jxo 

x €[a,b], n=l,2,.... (2.2) 

Remark 2.2. As we show below the introduced families of functions are closely 
related to the one-dimensional Schrodinger equations of the form u" — qu = 
Xu where q is a complex-valued continuous function. Slightly more general 
families of functions can be studied in relation to Sturm-Liouville equations 
of the form (py')' + qy = Xry. Their definition based on a corresponding 
recursive integration procedure is given in [42], [45], [37] . 

We introduce the infinite system of functions {</?fc} fe=0 defined as follows 

ff(x)XW{x), fcodd, 
<Pk(x) = < ,, , ?fB , . , (2.3) 

I f(x)X w (x), k even. 

The system (12.31) is closely related to the notion of the L-basis intro- 
duced and studied in [23] . Here the letter L corresponds to a linear ordinary 
differential operator. 

Together with the system of functions (|2.3j) we define the functions 
{V'fcjfcJLo usm g the "second half of the recursive integrals (|2.ip and (|2.2[) , 

,/ x , k odd, 
'■*('•) < X WL ( 2 - 4 ) 

,, s , k even. 

The following result obtained in |41) (for additional details and simpler proof 
see [32] and [35J) establishes the relation of the system of functions {fk}kL 
and {V'felfeLo t° the Sturm-Liouville equation. 

Theorem 2.3 (41 ). Let q be a continuous complex valued function of an 
independent real variable x £ [a, b] and X be an arbitrary complex number. 
Suppose there exists a solution f of the equation 

f" -qf = (2.5) 
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on (a, 6) such that f G C 2 (a,b) n C 1 ^, 6] and f(x) ^ for any x £ [a,b]. 
Then the general solution u £ C 2 (a, b) n C 1 ^, b] of the equation 

u" -gu = Xu (2.6) 

on (a, b) has the form 

U = C\U\ + C 2 U2 

where c\ and c 2 are arbitrary complex constants, 

oo ^ k oo ^ fe 

u i=E(^)j^ and U2 =J2^J-^^ k+ i (2.7) 

and both series converge uniformly on [a, b] together with the series of the 
first derivatives which have the form 



i 



f' + ^T^)\ ( -J<P2k+2kip 2 k-i ) and 



fe=l 



A fc //' 



^ (2fc + l)! V f 



<p 2k+l + (2k + l)ip 2k • (2 



T/ie series of the second derivatives converge uniformly on any segment 
[01,61] C (a, &). 

The representation (|2.7j) offers the linearly independent solutions of ()2.6p 
in the form of spectral parameter power series (SPPS). The possibility to 
represent solutions of the Sturm-Liouville equation in such form is by no 
means a novelty, though it is not a widely used tool (in fact, besides the work 
reviewed below and in [37] we are able to mention only [U Sect. 10], [33] and 
the recent paper [10]) and to our best knowledge for the first time it was 
applied for solving spectral problems in [35] . The reason of this underuse of 
the SPPS lies in the form in which the expansion coefficients were sought. 
Indeed, in previous works the calculation of coefficients was proposed in terms 
of successive integrals with the kernels in the form of iterated Green functions 
(see [U Sect. 10]). This makes any computation based on such representation 
difficult, less practical and even proofs of the most basic results like, e.g., the 
uniform convergence of the spectral parameter power series for any value of 
A 6 C (established in Theorem I2.3|) are not an easy task. For example, in [U 
p. 16] the parameter A is assumed to be small and no proof of convergence is 
given. 

The way of how the expansion coefficients in (I2.7P are calculated ac- 
cording to (|2.ip . (|2.2p is relatively simple and straightforward, this is why 
the estimation of the rate of convergence of the series (J2.7I) presents no dif- 
ficulty, see [IS]. Moreover, in [7J a discrete analogue of Theorem 12.31 was 
established and the discrete analogues of the series (|2.7j) resulted to be finite 
sums. 

Another crucial feature of the introduced representation of the expan- 
sion coefficients in (|2.7p consists in the fact that not only these coefficients 
(denoted by ipk in (|2.3p ) are required for solving different spectral problems 
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related to the Sturm-Liouville equation. Indeed, the functions J£( 2 *+ 1 ) and 
X^ 2k \ k = 0, 1, 2, . . . do not participate explicitly in the representation (|2.7p . 
Nevertheless, together with the functions tfk they appear in the representa- 
tion (12. 8|) of the derivatives of the solutions and therefore also in characteristic 
equations corresponding to the spectral problems. 

In the present work we also overview another approach developed in 
[43] . [8], [9] and [46] where the formal powers (|2.ip and (|2.2[) were consid- 
ered as infinite families of functions intimately related to the corresponding 
Sturm-Liouville operator. As we show this leads to a deeper understanding of 
the transmutation operators [3] , |10j also known as transformation operators 
[H] , [H] . Indeed, the functions ipt (x) result to be the images of the powers x k 
under the action of a corresponding transmutation operator [9] . This makes 
it possible to apply the transmutation operator even when the operator itself 
is unknown (and this is the usual situation - very few explicit examples are 
available) due to the fact that its action on every polynomial is known. This 
result was used in [5] and [S] to prove the completeness (Runge-type approx- 
imation theorems) for families of solutions of two-dimensional Schrodinger 
and Dirac equations with variable complex-valued coefficients. 

Remark 2.4. It is easy to see that by definition the solutions u\ and u 2 from 
(|2.7|) satisfy the following initial conditions 

ui(x ) = f(x ), u[(x ) = f(x ), 

u 2 (x ) = 0, u^o) = l/f(x Q ). 

Remark 2.5. It is worth mentioning that in the regular case the existence and 
construction of the required / presents no difficulty. Let q be real valued and 
continuous on [a, b]. Then (|2.5|) possesses two linearly independent regular 
solutions vi and vi whose zeros alternate. Thus one may choose / = v\ +iv2- 
Moreover, for the construction of v\ and v 2 in fact the same SPPS method 
may be used [35] . 



Theorem 12.31 together with the results on the completeness of Sturm- 
Liouville eigenfunctions and generalized eigenfunctions |51] implies the valid- 
ity of the following two statements. For their detailed proofs we refer to [33] 
and [44] respectively. 

Theorem 2.6 (43). Let (a, b) be a finite interval and f £ C 2 (a, b) n C l [a, b] 
be a complex valued function such that f(x) ^ for any x € [a, b]. 

If Xo = a (or xq — b) then each of the four systems of functions 

{* {2n) }Zo> {X (2n+1) }Zo> {X (2n) }Zo> {X (2n+1) }Zo « ™™ P lete l nL 2 (a,b). 

If Xq is an arbitrary point of the interval (a, b) then each of the fol- 
lowing two combined systems of functions {X' 2n H _ U |X' 2n+1 H and 
{X( 2n+1 )}~ =Q U {* (2n) } r 7 =0 is complete in L 2 (a, b). 

Theorem 2.7 ([44 ). Let f satisfy the conditions of the preceding theorem and 
{^Pk}kL be the system of functions defined by (|2.3[) with xo being an arbitrary 
point of the interval [a, b] . Then for any complex valued continuous, piecewise 
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continuously differentiable function h defined on [a, b] and for any e > there 
exists such N G N and such complex numbers ctk, k = 0, 1, . . . , N that 



max 

x£[a,b] 



N 

h{x) - y^afc^fc 

fc=0 



< e. 



3. Dispersion relations for spectral problems and approximate 
solutions 

The SPPS representation (|2.7p for solutions of the Sturm-Liouville equa- 
tion (|2.6j) is very convenient for writing down the dispersion (characteristic) 
relations in an analytical form. This fact was used in [13], [37], [35], |45) . 
[47] for approximating solutions of different eigenvalue problems. Here in 
order to explain this we consider two examples: the Sturm-Liouville prob- 
lem and the quantum-mechanical eigenvalue problem. As the performance of 
the SPPS method in application to classical (regular and singular) Sturm- 
Liouville problems was studied in detail in |45) here we consider the Sturm- 
Liouville problems with boundary conditions which depend on the spec- 
tral parameter A. This situation occurs in many physical models (see, e.g., 
0[Hl[T5J[2ni[17JEl] and references therein) and is considerably more difficult 
from the computational point of view. Moreover, as we show in this section 
the SPPS method is applicable to models admitting complex eigenvalues - 
an important advantage in comparison with the best purely numerical tech- 
niques all of them being based on the shooting method. 

Consider the equation u" — qu — \u together with the boundary condi- 
tions 

u(a) cos a + u (a) sin a = 0, (3-1) 

(3 lU (b) ~ &«'(&) = <f>(\)(P[u(b) - PW(b)), (3-2) 

where a is an arbitrary complex number, is a complex-valued function of 
the variable A and fix, /?2, (3[, f3' 2 are complex numbers. For some special 
forms of the function <j> such as <f>{\) = A or <fi(\) — A 2 + ciA + C2, results 
were obtained [19) . |64] concerning the regularity of the problem; we will 
not dwell upon the details. Notice that the SPPS approach is applicable as 
well to a more general Sturm-Liouville equation (pu')' + qu = Xru. For the 
corresponding details we refer to [37] and [45j . 

For simplicity, let us suppose that a = and hence the condition (|3.ip 
becomes u(a) — 0. Then choosing the initial integration point in (|2.ip and 
(J2.2I) as xq = a and taking into account Remark l2.4l we obtain that if an eigen- 
function exists it necessarily coincides with ui up to a multiplicative constant. 
In this case condition (|3.2p becomes equivalent to the equality [H], [37] 

(/w,(A)-/m(A))£^^«>( 6 >-^y^ x < 3 > >) = „, 

(3.3) 
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where ^2 (A) = 01$ — j3[ 2 4>W- This is the characteristic equation of the 
considered spectral problem. Calculation of eigenvalues given by (|3.3|) is espe- 
cially simple in the case of <fi being a polynomial of A. Precisely this particular 
situation was considered in all of the above-mentioned references concerning 
Sturm-Liouville problems with spectral parameter dependent boundary con- 
ditions. In any case the knowledge of an explicit characteristic equation (|3.3j) 
for the spectral problem makes possible its accurate and efficient solution. 
For this the infinite sums in p.3[) are truncated after a certain N £ N. 
The paper |45) contains several numerical tests corresponding to a variety of 
computationally difficult problems. All they reveal an excellent performance 
of the SPPS method. We do not review them here referring the interested 
reader to [45]. Instead we consider another interesting example from [37], a 
Sturm-Liouville problem admitting complex eigenvalues. 



Example 3.1. Consider the equation (|2.6p with q = on the interval (0,7r) 
with the boundary conditions u(0) = and u(ir) — —\ 2 u(n). The exact eigen- 
values of the problem are A„ = n together with the purely imaginary num- 
bers A± = ±i. Application of the SPPS method with N = 100 and 3000 inter- 
polating points (used for representing the integrands as splines) delivered the 
following results Ai = 1, A 2 = 4.0000000000007, A 3 = 9.00000000001, A 4 = 
15.99999999996, A 5 = 25.000000002, A 6 = 35.99999997, A 7 = 49.0000004, 
A 8 = 63.9999994, A 9 = 80.9996, A10 = 100.02 and A± = ±i. Thus, the com- 
plex eigenvalues are as easily and accurately detected by the SPPS method as 
the real eigenvalues. Note that for a better accuracy in calculation of higher 
eigenvalues of a Sturm-Liouville problem an additional simple shifting pro- 
cedure described in [JS] and based on the representation of solutions not as 
series in powers of A but in powers of (A — Ao) is helpful. We did not apply 
it here and hence the accuracy of the calculated value of A10 is considerably 
worse than the accuracy of the first calculated eigenvalues which in general 
can be improved by means of the mentioned shifting procedure. 

Figures 1-3 give us an idea about the stability of the computed eigen- 
values when N increases. In Fig. 1 we plot Ai and A2 computed with N = 
14, 16, ... , 120. Figure 2 shows A 3 computed with N = 24, 30, ... , 140 and 
Figure 3 shows A 4 computed with N = 40, 50, ... , 140 Similar plots can be 
done for calculated higher eigenvalues. In all cases the computed eigenvalues 
reveal a remarkable stability when N increases. 

An attractive feature of the SPPS method is the possibility to easily 
plot the characteristic relation. In Fig. 4 we show the absolute value of the 
expression from the left-hand side of (|3.3p as a function of the complex vari- 
able A for the considered example. Its zeros coincide with the eigenvalues 
of the problem. It is important to mention that such plot is obtained in a 
fraction of a second. This is due to the fact that once the required formal 
powers X^ are computed (and this takes several seconds) the calculation of 
the characteristic relation involves only simple algebraic operations. 
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FIGURE 1. The approximate eigenvalues Ai and A2 from Exam- 
ple [3TT] computed using different number N of formal powers. 
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FIGURE 2. The approximate values of A3 from Example l3.1l com- 
puted using different number N of formal powers. 
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FIGURE 3. The approximate values of A4 from Example l3.1l com- 
puted using different number N of formal powers. 
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FIGURE 4. The absolute value of the expression from the left- 
hand side of (|3.3p as a function of the complex variable A for the 
considered example. With the arrows we indicate the calculated 
complex eigenvalues A±. The other zeros of the graph correspond 
to the first five real eigenvalues of the problem. 
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Let us consider the one-dimensional Schrodinger equation 

Hu(x) = -u"{x) + Q(x)u(x) = \u{x), lei, (3.4) 

where 

{Qfl, x < 0, 

q(x), 0<x<h, (3.5) 

2 , x > h, 

ot\ and «2 are complex constants and q is a continuous complex-valued func- 
tion defined on the segment [0, h]. Thus, outside a finite segment the potential 
Q admits constant values, and at the end points of the segment the poten- 
tial may have discontinuities. We are looking for such values of the spectral 
parameter A G C for which the Schrodinger equation possesses a solution u 
belonging to the Sobolev space H 2 (R) which in the case of the potential of the 
form (|3.5|) means that we are looking for solutions exponentially decreasing 
at ±oo. This eigenvalue problem is one of the central in quantum mechanics 
for which H is a self-adjoint operator in L 2 (R) with the domain 7J 2 (R). It 
implies that Q is a real- valued function. In this case the operator H has a 
continuous spectrum [min {oi, 02} , +00) and a discrete spectrum located on 
the set 

min g(x), min {0:1,0:2} 
.xe[o,fc] 

Computation of energy levels of a quantum well described by the potential 
Q is a problem of physics of semiconductor nanostructures (see, e.g., |31|). 
Other important models which reduce to the spectral problem p. 41) arise in 
studying the electromagnetic and acoustic wave propagation in inhomoge- 
neous waveguides (see for instance [2], [16], [25], [E], [6], [56], [53]). 

A characteristic equation for this spectral problem in terms of spectral 
parameter power series was obtained in |13] (see also |37] ) where a simple nu- 
merical algorithm based on the approximation of the characteristic equation 
was implemented and compared to other known numerical techniques. Here 
we only give an example from |13j . 

The usual approach to numerical solution of the considered eigenvalue 
problem consists in applying the shooting method (see, e.g., [31]) which is 
known to be unstable, relatively slow and to the difference of the SPPS ap- 
proach does not offer any explicit equation for determining eigenvalues and 
eigenf unctions. In 30 another method based on approximation of the poten- 
tial by square wells was proposed. It is limited to the case of symmetric poten- 
tials. The approach based on the SPPS is completely different and does not 
require any shooting procedure, approximation of the potential or numerical 
differentiation. Derived from the exact characteristic equation its approxima- 
tion is considered, and in fact numerically the problem is reduced to finding 
zeros of a polynomial X)/c=o a k\^ m the interval [min<7(:z;),0), (/i 2 = —A). 

As an example, consider the potential Q defined by the expression 
Q(x) — — wsech x, x G (—00,00). It is not of a finite support, nevertheless 
its absolute value decreases rapidly when x — > ±00. The original problem is 
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approximated by a problem with a finite support potential Q defined by the 
equality 

{0, x < —a 

— v sech 2 x, —a < x < a 
0, x > a. 

An attractive feature of the potential Q is that its eigenvalues can be calcu- 
lated explicitly (see, e.g., [H]). In particular, for v — m(m+l) the eigenvalue 
A„ is given by the formula A„ = —(to — n) 2 , n = 0, 1, . . .. 

The results of application of the SPPS method for v = 12 are given in 
Table Q] in comparison with the exact values and the results from [30] . 



n 


Exact values 


Num. res. from [30] 


Num. res. 


using SPPS (N = 


-- 180) 





-9 


-9.094 


-8.999628656 


1 


-4 


-4.295 


-3.999998053 


2 


-1 


-0.885 


-0.999927816 



Table 1 . Approximations of \ n of the Hamiltonian 
H = -D 2 - 12 sech 2 x 



The results obtained by means of SPPS are considerably more accurate, 
and as was pointed out above the application of the SPPS method has much 
less restrictions. 



4. Transmutation operators 

We slightly modify here the definition given by Levitan [3S] adapting it to 
the purposes of the present work. Let E be a linear topological space and E\ 
its linear subspace (not necessarily closed). Let A and B be linear operators: 
Ex -> E. 

Definition 4.1. A linear invertible operator T defined on the whole E such 
that E\ is invariant under the action of T is called a transmutation operator 
for the pair of operators A and B if it fulfills the following two conditions. 

1. Both the operator T and its inverse T~ l are continuous in E; 

2. The following operator equality is valid 

AT = TB (4.1) 

or which is the same 

A = TBT- 1 . 

Very often in literature the transmutation operators are called the trans- 
formation operators. Here we keep ourselves to the original term coined 
by Delsarte and Lions [23 . Our main interest concerns the situation when 
A = —-7-5 + q(x), B = — ■4f, and q is a continuous complex- valued function. 
Hence for our purposes it will be sufficient to consider the functional space 
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E = C[a,b] with the topology of uniform convergence and its subspace E\ 
consisting of functions from C 2 [a, b] . One of the possibilities to introduce a 
transmutation operator on E was considered by Lions |50] and later on in 
other references (see, e.g., [51]), and consists in constructing a Volterra inte- 
gral operator corresponding to a midpoint of the segment of interest. As we 
begin with this transmutation operator it is convenient to consider a sym- 
metric segment [—a, a] and hence the functional space E — C[—a,a\. It is 
worth mentioning that other well known ways to construct the transmuta- 
tion operators (see, e.g., |49) . |63j ) imply imposing initial conditions on the 
functions and consequently lead to transmutation operators satisfying (|4.1|) 
only on subclasses of E\ . 

Thus, we consider the space E = C[—a, a] and an operator of transmu- 
tation for the defined above A and B can be realized in the form (see, e.g., 
[3S] and [5T]) of a Volterra integral operator 

Tu(x)=u{x) + K(x,t)u(i)dt (4.2) 

J — X 

where K(x,t) = H (^j^, ^j^) and H is the unique solution of the Goursat 
problem 

^^- = q(u + v)H(u,v), (4.3) 

H(u,0) = ^f q(s)ds, H(0,v)=0. (4.4) 

If the potential q is continuously differentiable, the kernel K itself is the 
solution of the Goursat problem 

' d 2 \ d 2 

Tr^-q{x) )K(x,t) = — K(x,t), 



dx 2 Hy J ) y ' ' dt 2 ~ 

1 f x 

K(x,x) = - / q(s)ds, K(x,—x)=0. 

2 Jo 

If the potential q is n times continuously differentiable, the kernel K(x,t) 
is n + 1 times continuously differentiable with respect to both independent 
variables (see [51]). 

An important property of this transmutation operator consists in the 
way how it maps solutions of the equation 

v" + lo 2 v = (4.5) 

into solutions of the equation 

u" — q(x)u + to u = (4-6) 

where w is a complex number. Denote by eo(iuj,x) the solution of (|4.6I) sat- 
isfying the initial conditions 

eo(iw,0) = 1 and e o (iw,0) = iio. 

The subindex "0" indicates that the initial conditions correspond to the point 
x = and the letter "e" reminds us that the initial values coincide with the 



14 Vladislav V. Kravchenko and Sergii M. Torba 

The transmutation operator (|4.2[) maps e lux into eo(iui,x), 

e (iu,x)=T[e iux ] (4.7) 

(see [EEJ Theorem 1.2.1]). 

Following [5T] we introduce the following notations 

K c {x, t-h) = h + K(x, t) + K(x, -t) + h I {K(x, £) - K(x, -£,)} d Z 

where h is a complex number, and 

K s {x, t; oo) = K(x, t) - K(x, -t). 

Theorem 4.2 ([51]). Solutions c(w,x;h) and s(u>,x; oo) of equation (|4.6|) sat- 
isfying the initial conditions 

c(u,0',h) = l, c' x (uj 1 0;h) = h (4.8) 

s(w,0;oo) = 0, 4(w,0;oo) = l (4.9) 

can &e represented in the form 

r 

c(uj,x;h) = cos u>x + / K c (x,t;h) cos uitdt (4-10) 

Jo 

and 

sinux f x „ , „ sinwi 



s(w,x;oo) = h / K s (x,t;oo) dt (4-11) 

Denote by 

T c u(x)=u(x)+ K c (x,t;h)u(t)dt (4.12) 

Jo 

and 

T s u(x)=u(x)+ K s (x,t;oo)u(t)dt (4.13) 



the corresponding integral operators. As was pointed out in |3], they are not 
transmutations on the whole subspace E\ , they even do not map all solutions 
of (|4.5j) into solutions of (|4.6|) . For example, as we show below 

-^ +qix) ) Ts[l] ^ Ts l~-^ {1) } =0 
when q is constant. 

.2 

Example 4.3. Transmutation operator for operators A := -pj + c, c is a 
constant, and i? := -^j. According to (|4.3p and (|4.4[) . finding the kernel 
of transmutation operator is equivalent to finding the function H(s,t) = 
X(s + 1, s — t) satisfying the Goursat problem 

The solution of this problem is given by [55J (4.85)] 

H(s,t) = --J J (2Vct(a-0)d? = -- ^ "> 
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where Jo and J\ are Bessel functions of the first kind, and the formula is 
valid even if the radicand is negative. Hence, 



R{ fx + y^ = _l_VME^WVME^l. (4.14) 

v ,y ' V 2 '2 J 2 x-y y ' 

From (J4.14I) we get the 'sine' kernel 



K s (x,t;co) = 5-^2 

and can check the above statement about the operator T s , 



r „ , f* tJc(x 2 - t 2 )JAJdx 2 - t 2 )) , ^ 

T 8 [l](x) - 1 - / -^ , V 2 V IJ -dt = J (a:v^, 

Jo a; - * 

d^ + c ) Ts[1] = x * - 

For the rest of this section suppose that / is a solution of (I2.5[) fulfilling 
the condition of Theorem I2.3I on a finite segment [—a, a]. We normalize / 
in such a way that /(0) = 1 and let /'(0) = h where h is some complex 
constant. Define the system of functions {fk}k=a by ^ ms function / with the 
use of (|2.1j) . (12. 2|) and (12. 3|) . The system of functions {fk}kL is related to 
the transmutation operators T c (with the same parameter h in the kernel) 
and T s in a way that it is the union of functions which are the result of acting 
of operator T s on the odd powers of independent variable and of operator T c 
on the even powers of independent variable. The following theorem holds, see 
[§] for the details of the proof. 

Theorem 4.4 (|9J). Let q be a continuous complex valued function of an in- 
dependent real variable x £ [—a, a], and f be a particular solution of (|2.5p 
such that f £ C 2 {—a, a), f ^ on [—a, a] and normalized as /(0) = 1. Let 
fk, k £ No := NU {0} are functions defined by (|2.3[) . Then the following 
equalities are valid 

fk = Tc[x "} when k £ No is even or equal to zero 

and 

<Pk = T s [x '} when k £ N is odd. 

As for the transmutation operator T, it does not map all powers of the 
independent variable into the functions (fk- Instead, the following theorem 
holds. 

Theorem 4.5 ([9 ). Under the conditions of Theorem \4--4\ the following equal- 
ities are valid 

(fk = T[x ] when k is odd (4-15) 

and 

ipk — -, fk+i = ^I 2 -"] when k £ No is even or equal to zero (4.16) 

rZ "T" ± 

where by h we denote /'(0) £ C. 
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Taking into account the first of former relations the second can be writ- 
ten also as follows 



•~Pk=T 



x k + J^x k+1 



when fc S Nq is even or equal to zero. 



fc + 1 

Remark 4.6. Let / be the solution of (|2.5p satisfying the initial conditions 

/(0) = 1, and /'(0) = 0. (4.17) 

If it does not vanish on [—a, a] then from Theorem 14.51 we obtain that 
(fk = T[x k ] for any k £ No- In general, of course there is no guaranty that 
the solution satisfying ()4.17|) have no zero on [—a, a]. Hence the operator T 
transmutes the powers of x into ifk whose construction is based on the solu- 
tion / satisfying (I4.17|) only in some neighborhood of the origin. In the next 
section we show how to change the operator T so that the new operator map 
x k into (flk(x) on the whole segment [—a, a]. 



Note that in Theorem [13] the operator T docs not depend on the func- 
tion /, so the right-hand sides of the equalities (|4.15|) and (|4.16[) do not 
change with the change of /. Consider two non- vanishing solutions / and g 
of (|2.5j) normalized as /(0) = g(0) = 1 and let ip k and ip k be the functions 
obtained from / and g respectively by means of (|2.1j) . (|2.2p and (|2.3|) . The 
relation between ip k and cpi are given by the following proposition which may 
be easily deduced from equalities f|4.15[) and (|4.16[) . 

Proposition 4.7. The following equalities hold 

ip k — ip 9 k when k £ N is odd, 

and 

rk = 'fit "I — i ~fk+i when k £ No is even, 

where hf = /'(0) and h g = g'(0). 

5. A parametrized family of transmutation operators 

In [S] a parametrized family of operators T/j, h £ C was introduced, given 
by the integral expression 

T h u(x) = u(x) + / K(x,t;h)u(t)dt (5.1) 

J — X 

where 

K(x, t;h)='^+ K{x, t) + J / (K{x, s) - K(x, -a)) ds. (5.2) 

2 2 J t 

They are related to operators T s and T c (with the parameter h in the kernel 
of the latter operator) by 

T h = T c P e + T S P , (5.3) 

where P e f(x) = (f(x) + f(-x))/2 and P f(x) = (f(x) - f(-x))/2 are 
projectors onto even and odd functions, respectively. In this section we show 
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that the operators T^ are transmutations, summarize their properties and in 
Theorem 15.81 we show how they act on powers of x. 

Let us notice that K(a;, t;0) = K(x,t) and that the expression 

K(x, t; h) - K(x, -t; h) = K(x, t) - K(x, -t)- 



(K(x, s) - K(x, -s)) ds = K(x, t) - K(x, -t) 



does not depend on h. Thus, it is possible to compute K(x, t; h) for any h by 
a given K(a;,£;/ii) for some particular value h±. We formulate this result as 
the following statement. 

Theorem 5.1 (|9J). The integral kernels K.(x,t;h) and K(sc, t\h\) are related 
by the expression 

K(x,t;h) = ^-+K(x,t;h l ) + ^^ J 



(K(x, s; hi) — K.(x , — s; hi)) ds. 

(5.4) 



The operator T/j may be expressed in terms of another operator Th t 
and in particular, in terms of the operator T. The following proposition holds. 

Proposition 5.2. The operators T; u and T/j 2 are related by the expression 

hi - hi 



Tfc, u = TV 



u(x) 



u(t) dt 



valid for any u G C[—a, a]. In particular, 



T,u = T 



u(x) + - J u(t) dt 



(5.5) 



(5.6) 



Proof. Using formulas (|5.1|) and (|5.2[) we obtain 

h r h r r 

T h u = Tu+~ u(t)dt+- u(t) K(x,s)dsdt- 
2 J-x 2 J_ x J t 



(t) / K(x,s)dsdt, 



h 

2 

and after changing the order of integration in the last two integrals, we have 



ThU = Tu - 



u(t) dt ■ 

h 

2 



K(x,s) / u(t)dtds 



h 
Tu + - 



u(t) dt+- J K(x, s) 



K(x,s) 



K(x,s) / u(t)dtds- 



u(t) dt 



u(t) dt ds = 



ds- 
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Tu+- I u(t)dt+- / K(x,s) J u{t)dtds = T 



h 

u(x) + - 



u{t) dt 



Since J_ j. u(s) ds dt = for any function u € C[—a, a], we have from (|5.6[) 
that 



■In 



In 

2 



u(x) 



U (ty 



u{t) dt 



hi 



= t u(x) ^ 

u(s) ds I dt 



lii 



u(t) dt- 



T 



u(x)- 



u(t) dt 



T h „u. 



D 



Using f|4.8[) f|4.13[) and (|5.3p it is possible to check how the operators 
T/j act on solutions of (14. 5t . 



Proposition 5.3 (|46j). The operator Th maps a solution v of an equation 
v" + lu 2 v = ; where u> is a complex number, into a solution u of the equation 
u" — q(x)u + (jj 2 u = with the following correspondence of the initial values 

u(0)=v(0), u'{0)=v'(0) + hv(0). (5.7) 

Remark 5.4. Formulas (|5.7[) are valid for any function v £ ^[—a, a]. 

We know that the integral kernel of the transmutation operator T is 
related to the solution of the Goursat problem (|4.3j) - (|4.4j) . A similar result 
holds for the operators T\. 

Theorem 5.5 ([46 ). In order for the function K{x, t; h) to be the kernel of a 
transmutation operator acting as described in Proposition \5.3\ it is necessary 
and sufficient that H(u, v; h) :— K(u + v, u — v; h) be a solution of the Goursat 
problem 

d 2 H{u,v;h) 
dudv 

\ + \j 9(s)ds, H(fi,v;h) = ± 

If the potential q is continuously differentiable, the function K(x,t;h) itself 
must be the solution of the Goursat problem 

d 2 , \ „, , s d 2 



H(u,0;h) 



= q(u + v)H(u, V] h), 

u 

q(s)ds, H(0,v,h) 



dx 2 



q(x)\K{x,t-h) 



dt 2 



K(x,t; h), 



K{x,x;h) = - + - 



q(s) ds, K(x, —x; h) — 



(5.8) 



(5.9) 



Under some additional requirements on the potential q the operators 
T/> are transmutations in the sense of Definition 14.11 The following theorem 
generalizes the results obtained in j46j . 

Theorem 5.6. Suppose the potential q satisfies either of the following two 
conditions. 

• q € ^[—a^]; 
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q G C[—a, a] and there exists a particular complex-valued solution g of 

(|2.5[) non-vanishing on [—a, a] . 

Then the operator T/j given by (|5.1j) satisfies the equality 



■—+q{x))T h [u} = T h 



d 2 
dx 2 



(5.10) 



for any u G C 2 [—a,a\. 

Proof. In [JS] the theorem was proved for the case q € ^f— a, a] and for the 
case when the particular solution g from the statement satisfies conditions 
g(0) = 1 and g'(0) = h. 

We may normalize the particular solution g as g(0) = 1. Suppose that 
#'(0) = hi. We know already that (|5.10[) holds for the operator T/^. To finish 
the proof, we use (|5.5[) and obtain 



.^ +ff (x))T fc [«]-(-^ + ? (x)]T /(1 



-hi 



-hi 



«"(a:) - 
u"(x) 



h-hi d 2 

2 d^ 



hi 



u{x) + 
u(t) dt 



h - hi 



i(t) dt 



u"{t) dt 



dx 2 



(") 



a 



Remark 5.7. As was pointed out in Remark [231 in the regular case the non- 
vanishing solution g of (|2.5|l exists due to the alternation of zeroes of two 
linearly independent solutions. Of course, it would be interesting to prove 
that the operators T^ are transmutations in the general case of complex- 
valued potentials q G C[—a, a] without any additional assumption. 

Suppose now that a function / is a solution of (|2.5|) . non- vanishing on 
[—a, a] and normalized as /(0) = 1. Let h :— /'(0) be some complex constant. 
Define as before the system of functions {ifik}kLo by this function / and by 
(|2.3p . The following theorem states that the operator T^ transmutes powers 
of x into the functions </?&. 

Theorem 5.8 ({9 ). Let q be a continuous complex valued function of an inde- 
pendent real variable x G [—a, a], and f be a particular solution of (|2.5|) such 
that f G C 2 {—a,a) together with 1// are bounded on [—a, a] and normalized 
as /(0) = 1, and let h := /'(0), where h is a complex number. Then the 
operator (|5.1|) with the kernel defined by (|5.2p transforms x k into ifik(x) for 
any k G No- 

Thus, we clarified that the system of functions {ifk} may be obtained 
as the result of the Volterra integral operator acting on powers of the in- 
dependent variable. As was mentioned before, this offers an algorithm for 
transmuting functions in the situation when K(ic, t',h) is unknown. More- 
over, properties of the Volterra integral operator such as boundedness and 
bounded invertibility in many functional spaces gives us a tool to prove the 
completeness of the system of function {tfk} in various situations. 
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Example 5.9. Consider a function k(x,t) = 2 f~+r\ a ter, in Example 16.81 it 
is explained how it can be obtained). We have 

{d 2 x - d 2 t )k{ X ,t) = - - 1)3 = j-^ ■ 2{x+l y 
k(x,—x) = 2(x+l) = ~2 an< ^ k(x,x) = 2^+T) = ~2 + 2 Jo [5+ip ^ S > thus the 

function k(x,t) satisfies the Goursat problem (I5.8[) - (|5.9I) with 
g(x) = 2/(a; + l) 2 and h — —1 and by Theorem 15.51 is the kernel of the 
transmutation operator T_i. 

Consider the function / = T_i[l] = —^j as a solution of (|2.5[) such that 
/(0) = 1 and /'(0) = h = —1, nonvanishing on any [—a, a] C (—1,1). The 
first 3 functions (fk are given by 

1 a; 3 + 3a; 2 + 3a; _ 2a; 3 + 3a; 2 

^ = /= ^TI' Vl= 3(3 + 1) ' ^" 3(a; + l) " 

It can be easily checked that 

r (t - l)t , a; 3 + 3a; 2 + 3a; 
T_ix = x + — — — - dt = — - — — r — = <p u 
J_ x 2(x + 1) 3(ar+l) 

2 /■* (t - l)t 2 ^ _ 2a; 3 + 3a; 2 



T_ix 2 = x l + / ■+ ^— di = — = ip 2 . 

^J_ x 2(x+l) 3(x+l) ^ 2 

We can calculate the kernel if of the original operator T by (|5.4p , it is given 

by 

T .. ± . 2x + 2t + x 2 - t 2 

K[X,t) — ; 

and we can check that T[x] = (ft and T[l] = - + ^f x ^ +3 — <^o + <Pi in 
accordance with Theorem 



6. Transmutation operators and Darboux transformed 
equations 

To construct the system of functions {fk}kLo we use the half of the functions 
{X( k \ XW}™ =Q . What about the second half? Note that starting with the 
function 1// wc obtain the same system of functions {X^ k \X^ k ^ with 
the only change that X\ becomes X^,\ and Xi ' becomes X^J f . Hence the 

"second half of the functions {JfW, X^}^ =Q from dH3J) is used. The func- 
tion 1// is continuous complex-valued and non- vanishing, and is a solution 

of the equation u" — qiu = 0, where q% = 2 I 4- ) — q. The last equation is 
known as the Darboux transformation of the original equation. The Darboux 
transformation is closely related to the factorization of the Schrodinger equa- 
tion, and nowadays it is used in dozens of works, see, e.g., pij 29, 52, 58] in 
connection with solitons and integrable systems, e.g., [TJ [321 [351 [57] and the 
review [59] of applications to quantum mechanics. 
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For the convenience denote the potential of the original equation by 
q\ and the corresponding Sturm-Liouville operator by A\ := ■4-^ — qi(x). 
Suppose a solution / of the equation A±f = is given such that f(x) ^ 
0, x € [—a 1 a] 1 it is normalized as /(0) = 1 and h := f'(0) is some complex 
number. Denote the Darboux-transformed operator by A2 := 4-% — q2(x), 

where q 2 (x) = 2(jffi) - qi(x). 

From the previous section we know that there exists a transmutation 
operator Ti ; /, for the original equation (|2.6[) with the potential q\ and such 
that 

T 1 . h x k = i p k , fee N . (6.1) 

The subindex "1" in the notation Ti-.h indicates that the transmutation op- 
erator corresponds to A\. 

Similarly, there exists a transmutation operator T-i;—h f° r the Darboux- 
transformed operator Ai such that 

T 2 ,- h x k = Vfc, k e No, (6.2) 

where the family of functions {V'fclfeLo i s defined by (|2.4|) . 

It is interesting to obtain some relations between the operators Ti-h 
and T2 ; _h and between their integral kernels Ki and K2. In this section we 
explain how to construct the integral kernel K2 by the known integral kernel 
Ki and show that the operators Ti ; h and T2-,-h satisfy certain commutation 
relations with the operator of differentiation. 

We remind some well known facts about the Darboux transformation. 
First, 1// is a solution of A 2 u = 0. Second, it is closely related to the 
factorization of Sturm-Liouville and one-dimensional Schrodinger operators. 
Namely, we have 

At = ^ - qi (x) = (d x + £) (d x -j) = jd x fd x j, (6.3) 

—,- q2 {x) = {d x -L)(d x + L)=fd^ 2 

Suppose that u is a solution of the equation A\u = uu for some ui G C. 
Then the function v — (d x — y) u = (f^xjju is a solution of the equation 
A2V = lov, and vice versa, given a solution v of ^v = wv, the function 
?i = (d x + y) v ~ (jdxf)v is a solution of Aiu = uiu. 

Suppose that the operator Ti := Ti ; /j which transmutes the operator 
A\ into the operator B — d 2 /dx 2 and the powers x k into the functions ifik is 
known in the sense that its kernel K.x(x,t;h) is given. As before h = /'(0). 
Then the function 1// is the non- vanishing solution of the equation A 2 u = 
satisfying l//(0) = 1 and (l//)'(0) = —h. Hence we are looking for the 
operator T2 := T2;-h transmuting the operator A2 into the operator B and 
the powers x k into the functions ipk- 

Let us explain the idea for obtaining the operator T2. We want to find an 
operator transforming solutions of the equation Bu + lu 2 u — into solutions 
of the equation A211 + ui 2 u — 0, see the first diagram below. Starting with 



A 2 = 3-5 - q 2 (x) = (d x - J - ){d x + J - )= fd x -d x f ■ . (6.4) 
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a solution a of the equation (d 2 + uj 2 )a = 0, by application of Ti we get a 
solution of {Ax + uj 2 )u — 0, and the expression (fd x j)Ti<j is a solution of 

(A 2 +uj 2 )v = 0. But the operator [fd x j)Ti is unbounded and hence cannot 
coincide with the operator T2 . In order to find the required bounded operator 
we may consider the second copy of the equation (d 2 + ui 2 )u = 0, which is a 
result of the Darboux transformation applied to {d 2 +lu 2 )<j = with respect 
to the particular solution g = 1 and construct the operator T2 by making 
the second diagram commutative. In order to obtain a bounded operator T2, 
instead of using fd x h for the last step, we will use the inverse of jd x f, i.e. 
j(f*f(s)-ds + C). 



Ti 


+ d* x 


- qi+uj 2 


d 2 x +L0 2 


*■%>- 


?1 


+ u: 2 


T^ 






fa.) 


d x 


j"9x/ )(ff-+C) 




dl 


-3 


2 + w 2 


d 2 x + u 2 


T2 a2 


<12 


+ 0J 2 



That explains how to obtain the following theorem. 

Theorem 6.1 ( 46J). The operator T2, acting on solutions u of equations {d 2 
Ld 2 )u = 0, lo G C by the rule 

1 



T 2 [u]{x) = 



f{x> wo 

coincides with the transmutation operator T2;-h 



/(7 ? )T 1 [ U '](7 7 )d7 ? + U (0) 



(6.5) 



Now we show that the operator T2 can be written as a Volterra integral 
operator and, as a consequence, extended by continuity to a wider class of 
functions. To obtain simpler expression for the integral kernel 1 K.2(x,t; —h) 
we have to suppose that the original integral kernel Ki(x,t; h) is known in 
the larger domain than required by definition (15.11) . Namely, suppose that the 
function Ki(x, t; h) is known and is continuously differentiable in the domain 
n : —a < x < a, —a < t < a. We refer the reader to 23 f° r further details. 

Theorem 6.2 ([46 ). The operator Ti admits a representation as the Volterra 
integral operator 



T 2 [u](x) = u{x) 



K 2 (x,t;-h)u(t)dt, 



with the kernel 



K 2 (x,t;-h) = -j^-(J X d t K 1 (s,f,h)f(s)ds+^f(-t) 

Such representation is valid for any function u £ C 1 [—a, a). 



(6.6) 



(6.7) 



By Theorems 16.11 and 16.21 the Volterra operators T2 and T2 coincide on 
the set of finite linear combinations of solutions of the equations (d 2 +lu 2 )u = 
0, u) G C. Since this set is dense in C[—a,a], by continuity of T2 and T2 we 
obtain the following corollaries. 
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Corollary 6.3 ( 46]). The operator T 2 given by (|6.6p with the kernel (|6.7[) 
coincides with T 2 on C[—a,a\. 

Corollary 6.4 ([46 ). The operator T 2 given by (|6.5[) coincides with T2 on 
C^-a.a}. 

Operator A\ is the Darboux transformation of the operator A 2 with 
respect to the solution 1//, hence we may obtain another relation between 
the operators Ti and T2. 

Corollary 6.5 ( 46J ). For any function u G C l [—a,a] the equality 

Ti[u](x) = f(x)(J j^-T 2 [u r }(r,)dv + u(0)) (6.8) 

is valid. 

From the second commutative diagram at the beginning of this subsec- 
tion we may deduce some commutation relations between the operators Ti, 
T2 and d/dx. The proof immediately follows from (J6.5I) and (J6.8I) . 

Corollary 6.6 (46 ). The following operator equalities hold on C l [—a,a\: 

d x fT 2 = fT x d x (6.9) 

d x jT 1 = jT 2 d x . (6.10) 

In [H] the following notion of generalized derivatives was introduced. 
Consider a function g assuming that both / and g possess the derivatives 
of all orders up to the order n on the segment [—a, a]. Then in [—a, a] the 
following generalized derivatives are defined 

7o(g)(x) =g{x), 

7 k (9)(x)=(f(x)) ( - 1)k ~ 1 ( 7k - 1 (g))'(x) 

for k = 1, 2, . . . ,?i. 

Let a function u be defined by the equality 

.9 = jTiU, 

and assume that u £ C n [—a, a]. Note that below we do not necessarily require 
that the functions / and g be from C n [—a, a]. With the use of (J6.9I) and (|6.10[) 
we have 

ll(g) = f ■ (j T l U )' = ^ ■ J T ^ U ' = / T 2 U '> 

12(g) = jz ■ (/Tau')' = -^ • /Tiu" = jTiu". 
By induction we obtain the following corollary. 
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Corollary 6.7 ( 46]). Let u E C n [-a,a] and g = JT lU . Then 

j k (g) = /T 2 M (fe) ifk is odd, k<n, 

and 

lk{g) — -jTiU^ ' ifk is even, k < n. 

Example 6.8. We start with the operator Aq = d 2 /dx 2 . We have to pick up 
such a solution / of the equation Aq/ = that /'// ^ 0. This is in order 
to obtain an operator A\ ^ Aq as a result of the Darboux transformation 
of A . For such solution / consider, e.g., fo(x) = X + 1. Both f and 1/fo 
are bounded on any segment [—a, a] C (— 1; 1) and the Darboux transformed 
operator has the form A x = -^ - ^jjp ■ 

The transmutation operator T for Aq is obviously an identity operator 
and Ko(x,t;0) = 0. Since /q(0) = 1, we look for the parametrized operator 
To ; i. Its kernel is given by (|5.4[) : Ko(a;,i;l) = 1/2. From Theorem 16.21 we 
obtain the transmutation kernel for the operator A\ 

Ki(M; -1) = —^ ■ ^ = 4^T V (6-H) 

x + 1 2 2{x + 1) 



the kernel from Example [ 

To obtain a less trivial example consider again the operator A\ = -4-^ — 
,2^ and the function f\{x) = (x + l) 2 as a solution of A\f = 0. Since 
h = f[(0) = 2, we compute Ki(x,t; 2) from (|6.11|) using ([5 

3x 2 + 6x + 4 - 3i 2 + 2t 



Ki(x,i;2) 



4(ar + l) 



The Darboux transformation of the operator .Ai with respect to the solution 
/i is the operator A2 — -j^ — fa-m 2 and by Theorem [621 the transmutation 
operator T2 ; _2 for A2 is given by the Volterra integral operator (|5.1j) with 
the kernel 



K2(M; - 2) =-(^w(£ 



2(.< 


!i+l , 


> + lf 


'■ds + 


(1-i) 2 


)- 




(3t- 


-l)(x 


+ 1) 2 


-3(i- 


l) 2 (t+l) 



4(x + l) 2 

This procedure may be continued iteratively. Consider the operators 

_ d 2 n(n+l) 
" := dx" 2 ~ (x + 1) 2 ' 

The function f n (x) — (x + l) n+1 is a solution of the equation A n f = 0. The 
Darboux transformation of the operator A n with respect to the solution /„ 
is the operator 

d 2 n (f^x)\ 2 , n(n + l) _ d 2 (n + l)(n + 2) 



dx 2 \ 



dx 2 \f n (x)J (x + 1) 2 dx 2 (x+iy 
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i.e., exactly the operator A n+ \. If we know K„(a>, t; —n) for the operator A n , 
by (|5.4[) we compute the kernel K n (x, i; n + 1) corresponding to the solution 
fn(%) and by Theorem 16.21 we may calculate the kernel K„ + i(x,t; — n — 1). 
Careful analysis shows that we have to integrate only polynomials in all 
integrals involved, so the described procedure can be performed up to any 
fixed n. 

Example 6.9. Consider the Schrodinger equation 

u" + 2sech 2 (ir)u = u. (6.12) 

This equation appears in soliton theory and as an example of a reflectionless 
potential in the one-dimensional quantum scattering theory (see, e.g. [48]). 
Equation (I6.12|) can be obtained as a result of the Darboux transformation 
of the equation u" = u with respect to the solution fix) — cosh a;. The 
transmutation operator for the operator A\ — <9 2 — 1 was calculated in [9, 
Example 3], Its kernel is given by the expression 



i Vx^¥h(Vx^¥) 

Kl(M ' 0) = 2 — ' 

where I\ is the modified Bessel function of the first kind. Hence from Theorem 
Owe obtain the transmutation kernel for the operator Ai = <9 2 +2 sech 2 x — 1 



K2( "' i;0) = 2^h^) J_ t { s-t + i^tr ) coshsds - 

7. Transmutation operator for the one-dimensional Dirac 
equation with a Lorentz scalar potential 

One-dimensional Dirac equations with Lorentz scalar potentials are widely 
studied (see, for example, [HI EH 121 121 EH [33 [Ml [Ml HH IE2] and [5S] for 
intertwining techniques for them). 

According to [54j the Dirac equation in one space dimension with a 
Lorentz scalar potential can be written as 

(d x +m + S{x))^> 1 =E^ 2 , (7.1) 

(-d x +m + S(x))^ 2 =E^ 1 , (7.2) 

where m (m > 0) is the mass and S(x) is a Lorentz scalar. Denote r\ — m + S 
and write the system (|7.ip . (17.21) in a matrix form as 

d x + v \ (*i\ _ E ( 1\ f^i 



o d m -nj ys>2) \-i oj V*2, 

In order to apply the results on the transmutation operators and factoriza- 
tions (|6.3p . (I6.4[) we consider a function / such that 

—— = -T}= -m-S{x). 
f(x) 

We can take f(x) = exp (— f* (m + S (s)) ds) , then /(0) = 1 and / does 
not vanish. Suppose the operators Ti and T2 are transmutations for the 
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operators A\ — (d x + 4-) (d x — 4-) and A 2 = (p x — 4-) (9 X + 4-) respectively 
(corresponding to functions / and 1// in the sense of Proposition 15. 3|) . As 
was shown in [46 with the use of commutation relations (|6.9p and (|6.10p , the 
operator 

'Ti 

,0 T 2/ 

transmutes any solution I of the system 

u[ = Eu 2 (7.3) 

u' 2 = -Em (7.4) 

into the solution ( ) of the system (|7.ip . (l7.2l) with the initial conditions 

Vw 2 y 

^1(0) = Ui(0), ^2(0) = "2(0). And vice versa if I T I is a solution of 
the system (|7.1[) , (|7.2[) . then the operator I 1 1 ) transmutes it into the 

solution ( Ul J of JI3]), CO} such that ui(0) = *i(0), u 2 (0) = * 2 (0). 

Consider two systems of functions {ipk}tLo an< ^ {V'fej^Lo constructed 
from the function / by (|2.3[) and (|2.4|l . The general solution of the system 
(|73 ]l .([7!3 ]l is given by 

Ml = C1W1 + C 2 f2 
M 2 = C*2Wl - C\V 2 , 

where C\ and C2 are arbitrary constants and 



v- (-WE' 
M x ) = cos Ex = 2_^ — (9 , v ■■- 



fc=0 



2fe 

(2*)! 



00 f 1UP2H1 



fe=0 



(2ft + 1)! 



From (|6.1j) and (|6.2p we see that the general solution of the one-dimensional 
Dirac system (|7.1j) . (|7.2[) has the form 

T „ ^ {-l) k E 2k „ ~ (_l)*£**+i 

$1 = Cl Z. (2fc)! ^ 2fe + C2 f. (2jfc + 1)! ^ +1 ' 

K/ — U hj — U 

Remark 7.1. It is possible to consider the two- or three-dimensional Dirac sys- 
tem and to construct the transmutation operator for it under some conditions 
on the potential. But the techniques involved, such as bicomplex numbers, 
pseudoanalytic function theory, Vekua equation and formal powers go well 
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beyond the scope of the present article. We refer interested readers to the 
recent paper [8]. 
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